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It is shown that behavior of the nuclear symmetry energy is the key quantity in the stability 
consideration in neutron star matter. The symmetry energy controls the position of crust-core 
transition and also may lead to new effects in the inner core of neutron star. 



INTRODUCTION 



In neutron star, the prevailing part of its interior is 
fulfilled with matter which is in the state called the beta 
equilibrium Already few meters below neutron star 
surface, at densities ^ lO^g/cm^ the degenerated elec- 
trons become relativistic and may easily penetrate the 
nuclei volume. The nuclcons, although being confined to 
nucleus, are subject to the beta equilibrium, e.i. they 
chemical potentials fulfill the constraint — Itp = Me- 
This leads to more and more neutron-rich nuclei and at 
the density pdri-p ~ lO^^g/cm'^ the neutron drip takes 
place. Above the neutron drip density one may consider 
the neutron star matter as a two-phase system in partial 
phase equilibrium. It is, indeed, because not all Gibbs 
conditions are fulfilled: P'^ = P""^ , = Me"'' . M^ = 
p'^^'^ but Mp Mp""' I where superscripts correspond re- 
spectively to the dripped and nuclear phase. At higher 
density, close to the saturation density the proton drip 
occurs as well Above this point p^ ~ /i™"^ f-iid 

one get a real two-phase system with all Gibbs condi- 
tions satisfied. The two phases have different properties 
like baryon density p and proton fraction x. Soon after 
the proton drip the differences in p and x vanishes and 
the star matter represents " npl matter" - a homogeneous 
nucleon-lepton liquid, (muons start to be produced at 
slightly higher densities). In the two-phase system, the 
nuclei form a lattice leading to a matter with solid state 
properties, unlike the homogeneous system being a fluid. 
In this way a neutron star has dense, liquid core covered 
by solid crust. The inclusion of finite-size effects, in some 
models leads to the presence of "funny phase" (rods, 
plates) in the transition region between the crust and 
core. In this short letter we limit ourselves to the bulk 
approximation and shall find the critical density where 
homogeneous npl matter splits into two phases not go- 
ing into details how the phase coexistence looks like. In 
this way we show that the crust-core transition is directly 
connected to the behavior of the nuclear symmetry en- 
ergy. The same analysis applied to higher density region 
reveal the same kind of instability and suggests repeated 
solidification in the central part of neutron star core. 



STABILITY CONDITIONS 

The beta reactions, taking place in neutron star, con- 
serve charge Q and baryon number B. Such statement 
means that any thermodynamical state in a volume V is 
uniquely described by setting the temperature T and all 
conserved quantity, here B and Q. We neglect the tem- 
perature effects which are relevant only for a very hot, 
proto-neutron star. Then the total energy U becomes a 
function of volume and conserved numbers U{V,B,Q). 
In order to consider stability of single phase one need 
to introduce intense (local) quantity u — U / B. The en- 
ergy per particle u becomes then a function of other lo- 
cal quantities, taken per baryon number v = V/B and 
q = Ql B. The first principle of thermodynamics takes 
the following form 

du = —P dv ~ p dq (1) 

where P is the pressure and p chemical potential of elec- 
tric charge. From beta equilibrium one may reads that 

p = Pe ^ Pl^L = Pn - Pp (2) 

The minus sign before p in (|^) comes from the defini- 
tion of Q = Np — Ne — which is negative for leptons 
Ip^ . The stability of any single phase, also called the 
intrinsic stability, is ensured by the convexity of u(v, q) 
Q. Thermodynamical identities allows to express this 
requirement in terms of following inequalities 1^ 

Usually, only the positive compressibility is examined, in 
particular, it is required for locally neutral matter that 

-m >o. ,4) 

However the second inequality in (^) is of the same im- 
portance. It concerns the stability of charge fluctuations 
and as it will be shown later it is connected to the positive 
value of the screening length in matter. Not all nuclear 
models ensure the charge fluctuations to be stable. As 



2 



was shown in the case of kaon condensation for wide class 
of models the system is not stable at any density [|j . In 
this work we would like to show that even such simple 
system like npl mater, without any exotic components, 
may represent a region of density where instability oc- 
curs. 

One may find another pair of inequalities which are 
equivalent to those in and, as it will appear later, are 
more convenient in further calculations: 
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(5) 



The intrinsic stability is determined by the details of 
nucleon-nucleon interactions. In order to show that, let's 
split the total energy per baryon into the nucleonic and 
leptonic part 



(6) 



The nucleonic contribution may be always expressed as 
a function of baryon number density n = B/V and the 
proton fraction x = Np/B. For leptons e^, the energy 
per volume, is completely determined by their chemical 
potential fi. Such decomposition is also true for the total 
pressure, so one may write 



u = (n, x) + e^{ii)/n 
P = P^{n,x)+P^{ji) 



(7) 
(8) 



where P = , henceforth, the partial derivatives 



of or P^ we indicate by subscripts n or x. One 
may show that Hn — IJ-p — {du^ / dx)n and then the beta 
equilibrium means that 



terms in the (^2|) refers to the nucleonic pressure and 
compressibility and they are positive from very funda- 
mental reasons. As we show later, the third term may 
lead to negative contribution and it comes from the lep- 
tons presence, it is the leptons that make the matter 
unstable. 

In order to find similar expression for the second in- 
equality in we make use of the expression for the 
total charge per baryon q ~ x ^ nj^/n , where the lepton 
number density rt^ depends on /i only (see Appendix). 
Differentiation of the equation by /i under constant n we 
obtain 



dq 
9/i 



dx\ 



(13) 



Here one may note that the second stability condition 
in (H) is directly connected to the stable screening of 
charged particles. The derivative appearing in ([T^), 
is proportional to the sum of squared inverse of screening 
lengths for leptons The first term in (|l3|), however, 
lacks this simple interpretation, it refers to hadronic in- 
teractions concealed by the beta equilibrium constraint. 
It is useful to call quantities like dq/dfi, as the "electri- 
cal capacitance" of matter. It measures the energetic cost 
of change in electric charge held in matter. For further 
discussion we introduce the compressibility and electric 
capacitance as 



K, 



dP 

dv 
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dn 

j ^ P,v 



i^q.li (14) 



(15) 



(9) 



Below we show that the equation (^) allows to express 
the stability conditions (||) or (^) in terms only of the 
nucleonic contribution to the total energy u. First, 
we take the compressibility from the second pair of in- 
equalities (0) 
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'dP 
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(10) 
(11) 



the leptons does not contribute when the derivative is 
taken under fixed. Fro m (p[ ) one get {dx/dn)f^ = 



,N 



/u^ and joining with (|ll[) we obtain 



dP 
dv 



,4 / 2 !^ 



-1-^1 (12) 



where for the right-hand side we applied the expression 
for the nucleonic pressure P^ 



n^u^ . The first two 



By use of the beta equilibrium relation (^ one may find 
that {dx/d^)n = ~^/u^x what allows to write the condi- 
tions (||) only in terms of nucleonic energy contribution 



,N 



K„ 



2 n + n'^u. 



H{ke + k^) 



,N 



> 0. 



> 



(16) 
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because also the second term in (O), coming from lep- 



tons, is determined by the relation ^ 



,N 



and nucle- 



onic pressure is = ri^u^^ . In the same way we may 
express the first pair of conditions (^ . In order to derive 
them it is useful to find identities between compressibil- 
ities and capacities. Then using standard theorem for 
implicit functions one may get the following relations 



XP -Xv = 
K„ 



Xva^n^ 
XP 

Xv 



(18) 
(19) 

(20) 
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with 



dn 



(21) 



Above relations perform analogy between the beta sta- 
ble matter and the one-component system at non-zero 
temperature. For instance, if one replace the chem- 
ical potential by the temperature, /i T, and the 
charge by the entropy g — > s in the Eg. ([l8|) , then re- 
cover the well known relations between specific heats: 
Cp~Cv ^ T{dPldV)T{dV/dTfp. Derivatives defined 
in ( pT| ) also may be expressed by means of (see Ap- 
pendix) so the conditions (|^) finally takes the form 



XP = 



2nu 



N 



2 N 



{2nL + n'^n u'^^) 



1 

,N „2 „,2 



> 



(22) 



> 



(23) 



For further analysis we choose the stability conditions 
expressed by the pair of equations (^,^). They are 
not only more readable than 23 ) but also more sensi- 
tive to the onset of instability. One may deduce it from 
Eg. ( p^ , p^ . Any stable system always keep 



and Xv < XP 



(24) 



N 



and thus the quantities if^ or Xv vanish first. 

As was mentioned above the nucleonic contribution u 
to the total energy, is a function of n and x. Isospin 
symmetry allows for the expansion in even powers of (1 — 
2x) but terms higher than quadratic are negligible, then 
the energy per baryon takes the form 



.N 



{n,x) = V{n)+Es{n){l - 2xf 



(25) 



where V{n) is the isoscalar potential and Eg (n) - the sym- 
metry energy, corresponding to the interactions isovector 
channel. Applying (||) to (|l|Jl^) we obtain 



= n^E'J{l-2x 



2(1 - 2xfEi^n^ 



Es 



Xv 



li{ke + k^) 



8E,in) 



V") +2n{E'g{l-2xf + V') 
> (26) 

> (27) 



Above equations shows explicitly the importance of sym- 
metry energy E^ in the stability considerations. The con- 
crete shape of the function Es{n) is not well known. We 
know only its value at saturation point, no = 0.16 fm~'^, 
where it takes about 30 MeV. Recent analysis on the 
isospin diffusion in heavy-ion collisions constrained sig- 
nificantly the slope of E'^{nQ) and the stiffness i?"(no) 
at saturation point ||^, ||] , however these results does not 



determine the high density behavior definitely. There are 
no experimental evidence about values of Es at very high 
density which is available in the central parts of a neutron 
star. In such extrapolations we must rely on the model 
calculations. For all of them the symmetry energy at sat- 
uration point have positive slope but at higher densities, 
they lead to different conclusions. For most the Eg is 
monotonically increasing function of n but some models 
lead to the Eg which saturates at higher densities or even 
bends down at some point and goes to zero ||l^, |l^, ^ . 
This kind of behavior is especially interesting as the last 
term in (^6|) may then take arbitrary large negative val- 
ues and lead to instability. From the other side, going 
to very low density, we encounter uncertainties as well. 
All models predict Eg decreasing to 0. However recent 
experiments show that symmetry energy saturates 
at the level about 10 MeV for very low densities. 



NUCLEAR MODELS 

In order to present the role played by the symmetry 
energy we apply a set nuclear models. To emphasize 
the Eg effects, the isoscalar part V{n) is kept the same, 
whereas the symmetry energy takes different forms. The 
isoscalar potential V{n), is that one taken from |^ which 
leads to the compressibility of nuclear matter equal to 
240 MeV at saturation point. For Eg we used two fam- 
ily of functions: four at low density and three for high 
density regime. At lower densities, we use shapes applied 
by Chen et.al. in there were numbered by a param- 
eter X = 1,0,-1,-2. Here we named them by a,b,c,d 
to avoid confusion with proton fraction x. The shapes of 
symmetry energy dependence at lower densities are pre- 
sented in Fig.0. At higher densities, much above np, we 
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FIG. 1: The different shapes of the symmetry energy at den- 
sities below saturation point. 

applied a "bent down" function. This type of symmetry 
energy with low values of Eg at high density is typical for 
realistic potentials pO|, pi] . Also in modern approaches 
like chiral dynamics ]12| and Skyrme effective forces |l^ 
small values of Eg were obtained. Here, for numerical 
simplicity, we introduced the simple polynomial function 
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which imitate results of works mentioned above 

E.H^Eio) <n-n,)in-n,) ^^^^ 
no (no - ni)(no - n2) 

where E^^^ = 30 MeV and n^ are the zeros of Eg (see 
Appendix). The shapes of these functions, named A, B, 
C are shown in the Fig.^. 
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FIG. 2: Three different shapes of the symmetry energy at 
densities above saturation point (solid lines). For comparison 
the results of realistic potentials (dotted lines), the higher: 
UV14UVII, and the lower: UV14TN1. 



LOW DENSITIES REGION 

The transition between liquid core and solid crust of 
the star is strictly connected to the breaking of the con- 
ditions (^, ^) . When at least one of them is broken the 
matter cannot be homogeneous anymore, it splits into 
two phases. The Fig.^ shows the compressibility under 
constant /x and its two contribution: " nuclear" if™"^ - the 
two first terms in ( p6| ) - and " beta" - the last term in 
( p6| ) which comes from the leptons presence. The " beta" 
contribution is always negative, hence there is always a 
competition between the positive " nuclear" compressibil- 
ity and the beta reactions which tends to destabilize the 
matter. At some critical point, ric, the compressibility 
changes its sign and below n,. the matter cannot exist 
as a single, neutral phase. The actual splitting into two 
phases does not occur exactly at n^, but slightly above 
TLc, because the system must find a state where the two 
charged phases may coexist. The correction is expected 
to be tiny, so the critical point for the vanishing of K^^ 
may be treated as a good estimation for the boundary 
of the liquid core in neutron star. The Table | shows 
the critical density at which if^ vanishes. It depends 
strongly on E^ but does not behaves monotonically with 
the Es- For symmetry energy taking both high and low 
values (a, d case) the Uc moves to higher density close 
to saturation point. The lowest values of Uc are achieved 
with intermediate Eg (models b and c). There is no sim- 
ple relation between values of Es and Uc because the first 
and the second derivatives of Eg are essential as well. 
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FIG. 3: The compressibility Tf^ (thick) and its contributions 
(thin lines). The dotted line corresponds to the energy per 
baryon for neutral matter u{n, 0). 



TABLE 1: The critical density for different models. 



model 


a b 


c d 


Wc, fm"'' 


0.119 0.092 


0.095 0.160 



HIGH DENSITIES REGION 

The stability of matter at densities much greater than 
no do not need to be obvious. For the symmetry energy 
increasing in the whole range of density the matter is sta- 
ble indeed. However, the chosen nuclear models with very 
low values of Eg lead again to the same kind of instability 
as occurs in the crust-core transition region. For all pre- 
sented models A,B,C there is a critical density ric where 
vanishes. The behavior of the compressibility for 
the model B is shown on the FigJJ. It is worth to notice 
that energy per baryon for neutral matter u{n) = u(n, 0) 
does not reveal any pathology, it is monotonically in- 
creasing and its convexity, u„„(n) = Kqjr? , remains pos- 
itive. If one look only on the energy per baryon behavior, 
one may overlook that the matter becomes unstable at 
some point. 
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FIG. 4: The compressibility (thick) and its contributions 
(thin lines). The dashed line corresponds to Kq and dotted 
hue energy per baryon for neutral matter, it(n) = ii(n, 0) 

Of course, the instability point has physical meaning 
only if it is attainable in a neutron star. The Tab.O 
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shows the neutron star properties (the central density 
tlcen of a star with maximal mass Mmax found by solving 
of TOV equation jl^. Actually, in the case A and B 

TABLE II: The critical density and neutron star properties 
for different models. 



model 


A 


B 


C 


Tlcen, fm"'' 


1.92 


1.32 


1.21 




1.64 


1.73 


1.84 


Uc, fm~^ 


0.74 


1.20 


1.43 



the Tic is smaller than the central density of maximal 
star. It means that for sufficiently massive stars their 
structure changes essentially in the central part of the 
star. The homogeneous phase again splits into two-phase 
system. One may suspect formation of funny phases or 
solidification of the central part of stellar core. 



The parameterization of Eg at high density regime: 



model 


A 


B 


C 


ni, fm"'' 


1.0 


1.5 


1.8 


n2, fm"^ 


2.3 


2.5 


10.0 



SUMMARY 



In this letter we would like to notice the role played by 
the nuclear energy symmetry Es in the stability of dense 
matter under beta equilibrium. It was shown that rele- 
vant quantity in such considerations is - the compress- 
ibility under constant chemical potential, rather than Kq 
- the compressibility under constant charge. The instabil- 
ity of matter leads to phase separations and corresponds 
to the transition from the liquid core to the solid crust. 
In this way one may get simple connection between Eg 
behavior at low densities and the size of the star crust, 
which may be estimated from pulsar glitching jist . 

The stability considerations were also carried out at 
very high density. It was shown that for nuclear mod- 
els with small values of E^ it is possible that instability 
occurs. This may lead to repeated solidification the in- 
ternal parts of a core. Such solid inner core would have 
important consequences on the rotation of a star and its 
magnetic properties. 



APPENDIX 

Here we present formulas mentioned in the regular 
text, the lepton number density 

and mixed derivatives appearing in and ( p^ 
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dn 
dn 



N 2 
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1 + U^x^'lV 
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The results of TOV equations are presented under as- 
sumption that single phase EOS is stable for all densi- 
ties up to the max density, also above ric. It only shows 
that ric is attainable for a given EOS. Above ric the EOS 
should be corrected and ricen may change but not much. 



